Effect of epitaxial strain on the spontaneous polarization of thin film ferroelectrics 
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We investigate the variation of the spontaneous ferroelectric polarization with epitaxial strain 
for BaTiOa, PbTiOa, LiNbOa, and BiFeOa using first principles calculations. We find that while 
the strain dependence of the polarization is very strong in the simple perovskite systems BaTi03 
and PbTiC>3 it is only weak in LiNbC>3 and BiFeOs. We show that this different behavior can be 
understood purely in terms of the piezoelectric and elastic constants of the unstrained bulk material, 
and we discuss several factors that determine the strain behavior of a certain material. 
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The possible application of ferroelectric materials in 
microelectronic devices has led to strong interest in the 
properties of thin film ferroelectrics One important 
question in this context is how epitaxial strain, which is 
incorporated in the ferroelectric material due to the lat- 
tice mismatch with the substrate, affects the ferroelec- 
tric characteristics of the thin film, ft has been demon- 
strated that epitaxial strain can have drastic effects, such 
as inducing ferroelectricity at room temperature in oth- 
erwise paraelectric SrTiC>3 0, or increasing the ferroelec- 
tric Curie temperature of BaTiOa by nearly 500°C and 
the remanent polarization by 250 % compared with the 
corresponding bulk values Based on these observa- 
tions it is often assumed that the strong sensitivity to 
epitaxial strain is a common feature of all ferroelectrics. 

Indeed, this assumption has been supported by first 
principles calculations for several simple ABO3 per- 
ovskite ferroelectrics 0, H, IE 0- Strain is introduced 
in these calculations by fixing the lattice constants corre- 
sponding to the lateral directions of the substrate while 
relaxing all remaining structural parameters, according 
to the minimum of the total energy under the epitax- 
ial constraint. This makes it possible to isolate the 
pure strain effect from other effects present in real thin 
film samples such as structural defects, chemical inhomo- 
geneities, and interface effects. For example, in Refs. |fj| 
and0 the change in phase stability due to epitaxial strain 
was investigated by first principles techniques for several 
simple perovskite systems including the prototype fer- 
roelectrics BaTiOa and PbTi03. It was shown that for 
these systems a series of consecutive phase transitions oc- 
cur which effectively rotate the polarization from "out- 
of-plane" for compressive epitaxial strain to "in-plane" 
for tensile epitaxial strain. Also, the magnitude of the 
spontaneous polarization within the different phases was 
shown to be strongly strain dependent. In contrast, we 
have recently shown that in BiFeOs, a multiferroic sys- 
tem where ferroelectricity coexists with magnetic order, 
the magnitude of the polarization barely changes even for 
relatively large strain values of ±3 % y|. 

In the present work we investigate the factors that lead 
to a strong strain dependence of the spontaneous polar- 
ization in some materials and to a relative inertness to 



epitaxial strain in other materials. To do this we com- 
pare the spontaneous polarization as a function of epi- 
taxial strain (in a fixed phase) for a variety of different 
perovskite-derived ferroelectrics with different structural 
symmetries and different mechanisms causing the ferro- 
electric displacements. We show that the sensitivity to 
epitaxial strain varies considerably for different materials 
and that strong strain dependence is not a universal fea- 
ture of all ferroelectrics. Furthermore, we show that the 
strain dependence of the polarization for experimentally 
relevant strain values in all systems can be understood in 
terms of the piezoelectric and elastic constants of the un- 
strained materials and we discuss the factors determining 
the sensitivity to epitaxial strain of a certain material. 

In general, the total change of the spontaneous polar- 
ization P to linear order in the strain components e$ is 
given by the improper piezoelectric tensor, c a i 



dP a 



(1) 



Here, a = 1,2,3 stands for the usual cartesian com- 
ponents x, y, and z, whereas i — 1,...,6 denotes the 
components of the strain tensor in Voigt notation (see 
e.g. Ref. E3)- 

For epitaxial strain, 64 = £5 = e§ = 
and £1 = £2 ^ 0. Furthermore, the ratio between in- 
plane and out-of-plane strain is given by the Poisson ra- 
tio 11 = —61/63. In the following we only consider the 
situation where the spontaneous polarization is directed 
perpendicular to the orientation of the substrate plane, 
and we use this as the z direction of our coordinate sys- 
tem. The substrate surface is then oriented parallel to 
the x-y plane. In this case, the change in polarization 
(to linear order in the strain) is given by: 



AP, 



(2c 3 i 



£33\ 

n ) 



ei = c c ff ei 



(2) 



Here we have further assumed that the symmetry is such 
that c 3 i = c 32 . 

We now investigate the validity of Eq.|21for experimen- 
tally relevant strain values by comparing directly calcu- 
lated values of the spontaneous polarization as a func- 
tion of epitaxial strain with values calculated using Eq.|2 
and the piezoelectric constants and Poisson ratios of the 
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corresponding unstrained structure. We also investigate 
how the effective piezoelectric constant c e ff, which de- 
scribes the change in polarization due to epitaxial strain, 
varies for different materials. 

In Ref. several possible explanations for the differ- 
ent strain behavior in BiFeOs compared to systems like 
BaTiOa/PbTiOs were suggested: (i) different structural 
symmetry: BiFe0 3 crystallizes in the rhombohedral i?3c 
structure with a doubled unit cell compared to the ideal 
perovskite structure and shows rotated oxygen octahe- 
dra in addition to the polar displacements (see Ref. ITlh . 
whereas the relevant phases [iij for BaTi0 3 /PbTi0 3 are 
simple tetragonally distorted perovskite structures (space 
group PAmm) 0, 0] ; (ii) different displacement mecha- 
nisms: in BaTiOs/PbTiOa the d° configuration of the 
Ti cation plays a crucial role for the ferroelectric insta- 
bility, whereas in the multiferroic system BiFc03 the fcr- 
roelectricity is solely caused by the Bi 6s lone electron 
pair [J^; (iii) a general high stability of the ferroelec- 
tric state in BiFcC>3, indicated by rather large ionic dis- 
placements and a high ferroelectric Curie t emp erature of 
1123 K (BaTi0 3 : 400 K, PbTi0 3 : 763 K) |ig. 

In order to systematically address the possible explana- 
tions (i) to (iii), we study the following systems: the pro- 
totypical ferroelectrics BaTi03 and PbTi0 3 ; the multi- 
ferroic system BiFeOs, both in its ground state i?3c struc- 
ture and in a hypothetical metastable PAmm structure, 
i.e. isostructural to BaTi0 3 /PbTi0 3 ; LiNb0 3 , which is 
isostructural to BiFeC>3, with large ionic displacements 
and a very high Curie temperature of 1480 K 0, but 
where, in contrast to the multiferroic system BiFeOs, the 
transition metal cation plays an active role in the ferro- 
electric instability, similar to the case of BaTi0 3 /PbTi0 3 

MM 

Values for the spontaneous polarization obtained in 
this work are calculated using the Berry phase approach 
for determining the electronic contribution to the polar- 
ization [l6l ITtI ] . The piezoelectric constants are obtained 
using the Berry phase approach and a finite difference 
method analogous to Ref. [l8l The Poisson ratios are es- 
timated from the change in the out-of-plane lattice con- 
stant when fixing the in-plane lattice constant and re- 
laxing all ionic positions. For all calculations we use the 
Vienna Ab-initio Simulation Package (VASP) 19], which 
employs the Projector Augmented Wave (PAW) method 
within density functional theory j^, • Except where 
otherwise noted we use the local spin-density approxima- 
tion (LSDA) , see Ref. H3- All calculational parameters 
are tested to result in good convergence of all quantities 
under consideration. 

The structural parameters calculated in this work cor- 
responding to the unstrained systems are listed in Tabled 
for LiNb0 3 and in Table ITl) for the tetragonal systems. 
All structural parameters are i n g ood agreement with 
previous calculations, see Refs. 1 14] and lla for LiNb03, 
Ref. for BaTi0 3 , and Ref. Hffor PbTi0 3 . For i?3c 



TABLE I: Structural parameters calculated for LiNbC>3. a 
and c are the lattice parameters of the corresponding hexag- 
onal unit cell. z,u, v, w are internal structural parameters 
defined in Ref. fl5t 



a [A] c [A] z 


u 


V 


10 


5.073 13.702 0.0337 


0.01257 


0.0426 


0.0178 


TABLE II: Structural parameters calculated for PAmm 
BaTi03, PbTiOa, and BiFeOs. a and c/a describe the tetrag- 
onal unit cell, Aj is the displacement of ion i along the c axis 
(in units of c). The notation is the same as used in Ref. 


a [A] c/a 


A B 


A 0/ 


Ao„ 


BaTi0 3 3.945 1.011 
PbTiOs 3.870 1.040 
BiFeOs 3.665 1.270 


0.012 
0.034 
0.064 


-0.019 
0.087 
0.170 


-0.013 
0.101 
0.203 



BiFeO s th e structural parameters are listed in Table II 
of Ref.TTll (U e ff = 0). All systems except PAmm BiFe0 3 



are fully relaxed using the LSDA. For PAram BiFe0 3 cal- 
culations are performed using the LSDA+[/ method [24| 
with a small [/ c ft = 2 eV in order to ensure the insulating 
character. The Poisson ratios and the spontaneous po- 
larization as a function of epitaxial strain were extracted 
from Ref. for BaTi0 3 , and from Ref. 6 for PbTi0 3 . All 
other quantities are calculated in this work. 

Since BiFe0 3 in the metastable PAmm structure has 
not been discussed in the literature before, we give a 
few more details obtained by our structural relaxation 
of this system. The energy minimum of the relaxed 
PAmm structure of BiFe0 3 is 0.26 eV/(formula unit) 
higher than for the ground state i?3c structure (calcu- 
lated with the same U e g). The relaxed PAmm structure 
is strongly distorted compared to the ideal perovskite 
structure, with a c/a ratio of 1.27 and a very large polar- 
ization of 151 /LtC/cm 2 . These values are similar to those 
found for the recently predicted isostructural ferroelectric 
BiGaOs M 

Before presenting our main results, we clarify two tech- 
nical points. First, since the polarization of a periodic 
solid is only well defined modulo eR/V, where e is the 
electronic charge, R is a lattice vector, and V is the unit 
cell volume, the calculation of the polarization for a cer- 
tain bulk structure leads to a lattice of values 0, 0] . 
The spontaneous polarization P can be obtained by mon- 
itoring the change in polarization AP of an arbitrary cho- 
sen "branch" of the polarization lattice, when the struc- 
ture is deformed from the original structure into a state 
with inverted polarization. It has been shown by King- 
Smith and Vanderbilt 0, 0] that this change is inde- 
pendent of the chosen "path" leading from the initial to 
the final state, as long as the system stays insulating for 
all intermediate states along the path. The spontaneous 
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Py=d/V X 



P 2 =d/V 2 <P l 



FIG. 1: Reason for the difference between proper and im- 
proper piezoelectric response: the polarization P2 is smaller 
than Pi although the dipole moment d is unchanged. 



TABLE III: Improper Piezoelectric constants C33, C31, and 
c e ff (in /xC/cm 2 ) as well as the Poisson ratio n for the various 
systems. References are given for values not calculated in this 
work. 





C33 


C31 


n 


Cefi 


BaTiOa 


670 


30 


0.65 [4] 


-971 


PbTi0 3 


586 


103 


0.58 [6] 


-804 


BiFeOs (i?3c) 


213 


50 


1.15 


-85 


LiNbOs 


124 


-75 


1.5 


-233 


BiFe0 3 (P4mm) 


203 


-15 


0.65 


-342 



150 



100 



u 



50 



BiFe0 3 (P4mm) - 


- PbTi0 3 


BiFeCL (R3c) 




LiNb0 3 




1 


^^aTi0 3 



-2 2 

e t [%] 



polarization is then given by P — AP/2 and is also in- 
dependent of the chosen branch. 

Second, it has been pointed out in the literature that 
the improper piezoelectric response depends on the spe- 
cific branch used to calculate this quantity [5fij . The 
distinction between proper and improper piezoelectric re- 
sponse is a consequence of the definition of the polariza- 
tion as the dipole moment per volume, which causes the 
polarization to vary as a function of strain even when the 
dipole moment does not change (see Fig. ^| 9, 26] . This 
pure volume effect is included in the improper but not in 
the proper piezoelectric response. However, the two are 
trivially related by the spontaneous polarization (in the 
case of C33 the proper is equal to the improper piezoelec- 
tric response, for C31 the proper piezoelectric response is 
given by the improper value plus the spontaneous polar- 
ization) [9II26I]. Here we report the improper piezoelectric 
response, which is used in Eqs. (JTJ and J2J, since they re- 
fer to the total change in spontaneous polarization. Wc 
point out that the improper piezoelectric constant corre- 
sponding to the spontaneous polarization is well defined 
and does not depend on the specific branch used in the 
evaluation of this quantity. 

The improper piezoelectric constants and Poisson ra- 
tios for the five systems are given in Table llHl Note that 
for the tetragonal systems the fourfold axis is used as 
the z direction of our coordinate system, whereas for the 
rhombohcdral systems the threefold axis is used. It is evi- 
dent that the effective piezoelectric constants c G ff describ- 
ing the effect of epitaxial strain (to linear order) on the 
spontaneous polarization are much larger for BaTiOa and 
PbTiOs than for the other three systems. This is mainly 
due to the larger C33 in BaTiOs/PbTiOs but also due to 



FIG. 2: (Color online) Spontaneous polarization P for sev- 
eral ferroelectrics as a function of epitaxial strain ei. Sym- 
bols represent directly calculated values, lines are calculated 
using Eq. and the effective piezoelectric constants from 
Table ETT1 Directly calculated values for BaTi0 3 and PbTi0 3 
are taken from Refs.0 and[(J respectively. 

the smaller Poisson ratio of the PAram structures [3fj|. 
The Poisson ratios are comparable for systems with the 
same structural symmetry. In PAmm BiFeOs the smaller 
Poisson ratio leads to a slightly increased c c g compared 
to R3c BiFeC>3 although the piezoelectric constant C33 is 
similar for both symmetries. Another observation is that, 
although C33 is smaller in Ric BiFeOs than in LiNbC>3, 
Coff is larger due to the different sign of C31 . 

Fig. |21 compares the directly calculated change in spon- 
taneous polarization caused by epitaxial strain with the 
corresponding change calculated using c c s from Table iTTll 
for all five systems. One can see that the linear approxi- 
mation of Eq. |2] works well for all systems up to strains of 
± 3 %. Even in the case of BaTi0 3 and PbTi0 3 , where 
we have used data from two different sources, the agree- 
ment between the polarization calculated using Eq.^and 
the directly calculated values is remarkably good. For 
the systems where all values are calculated in the present 
work the agreement is extremely good. 

The present study shows that the polarization response 
to epitaxial strain for experimentally relevant strain val- 
ues can be described in terms of the piezoelectric and 
elastic constants of the unstrained system. Two further 
observations can be made. First, the appearance of two 
different terms in Eq. [3 with potentially opposite signs 
makes it possible that a system can in principle have 
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rather large piezoelectric constants, but that due to par- 
tial cancellation c e g is rather small. Second, a large value 
of c 33 can eventually be suppressed by a large Poisson 
ratio. Otherwise, if no partial cancellation occurs, a sys- 
tem with a large piezoelectric response will in general also 
show a strong response of the spontaneous polarization 
to epitaxial strain. 

The question of whether epitaxial strain enhances the 
polarization in ferroelectric thin films therefore leads to 
the more general problem of what determines the piezo- 
electric (and elastic) response of a certain material. An 
overview of recent first principles research on piezoelec- 
tricity is given in Ref. |27j. Some insight can be gained 
by separating the piezoelectric response into a "clamped 
ion" part and an "internal strain" part. The first mea- 
sures the strain-induced change in polarization for fixed 
internal ionic coordinates whereas the latter gives the 
contribution resulting from the relaxation of the ions. It 
can be shown that the internal strain part of the piezo- 
electric tensor contains the inverse of the force constant 
matrix j^. Therefore, if the force constant matrix has 
rather small eigenvalues, i.e. low energy phonon modes, 
and these modes couple to a certain strain component, 
the piezoelectric response can be very large and is then 
dominated by the internal strain part. In a soft mode 
system at least one phonon energy goes to zero at the 
ferroelectric transition temperature and the piezoelectric 
response diverges . Far away from the transition tem- 
perature the energy of this soft phonon mode increases 
again. This can in principle explain the decrease of the 
(zero temperature) value of C33 in the series BaTiC>3, 
PbTiC>3, BiFc03, and LiNbC>3, i.e. for increasing Curie 
temperatures. For C31 both contributions to the piezo- 
electric constant are comparable and of opposite sign, 
leading to no clear trend in magnitude or sign for the 
different materials. 

In summary, regarding the three possible explanations 

(i) to (iii) given in Ref. |g for the different strain behavior 
in different systems, we can now conclude that the struc- 
tural symmetry (i) mainly affects the value of the Poisson 
ratio, whereas the "general stability" of the ferroelectric 
state (iii) determines the order of magnitude of C33. No 
clear influence of the different displacement mechanisms 

(ii) or of the "multiferroicity" of BiFe03 can be identified 
from our results, although the displacement mechanism 
can of course influence the stability of the ferroelectric 
state. It seems that in most cases a large C33 leads to a 
strong dependence of the polarization on epitaxial strain, 
and conversely a strong epitaxial strain dependence is 
usually connected with a large C33. Nevertheless, a small 
(large) Poisson ratio can enhance (suppress) this effect 
and the same is true for the sign and magnitude of C31 . It 
also becomes clear from our results that not every ferro- 
electric necessarily shows a strong sensitivity to epitaxial 
strain, and that the response for a certain material can 
be calculated from its piezoelectric tensor and Poisson ra- 



tio using Eq. [3] for all experimentally accessible epitaxial 
strains. 
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[29] BaTi03 is tetragonal in the temperature range 280 K < 
T < 400 K 0. 

[30] The piezoelectric constants for PbTiOa have also been 
calculated in Ref. using the generalized gradient ap- 
proximation (GGA) instead of the LSDA. It seems that 



the LSDA values are slightly larger than the GGA val- 
ues, but the difference could also (at least partly) be due 
to the different structural parameters used in the two 
calculations. 



